ac Josephson effect in asymmetric superconducting quantum point contacts 
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We investigate ac Josephson effects between two superconductors connected by a single-mode 
quantum point contact, where the gap amplitudes in the two superconductors are unequal. In these 
systems, it was found in previous studies on the dc effects that, besides the Andreev bound-states, 
the continuum states can also contribute to the current. Using the quasiclassical formulation, we 
calculate the current- voltage characteristics for general transmission D of the point contact. To 
emphasize bound versus continuum states, we examine in detail the low bias, ballistic (D = 1) 
limit. It is shown that in this limit the current- voltage characteristics can be determined from 
the current-phase relation, if we pay particular attention to the different behaviors of these states 
under the bias voltage. For unequal gap configurations, the continuum states give rise to non- 
zero sine components. We also demonstrate that in this limit the temperature dependence of the dc 
component follows tanh(A s /2T), where A s is the smaller gap, with the contribution coming entirely 
from the bound state. 

PACS numbers: 74.80.Fp, 74.50.+r 



I. INTRODUCTION 



When a thin normal layer separates two superconduc- 
tors, the superconducting coherence can spread across 
the normal region. The (dc and ac) Josephson effects 
are typical manifestations of such phenomenal The rise 
of mesoscopic physics leads to reconsideration of these 
effects in systems where the normal region consists of 
narrow channels with quantized transverse modes (quan- 
tum point contacts). Typically the superconductors are 
connected via quantum point contacts which can be con- 
strictions in semiconductor heterostructures^ or atomic 
contacts in break junctions^ One of the insights gained 
from these studies is the important role of the Andreev 
bound-states in the Josephson effects. 4 

When a particle (hole) incident from a normal metal 
into a superconductor (or vice versa), besides normal re- 
flections, the particle (hole) can be retro-reflected along 
its incident path and converted into a hole (particle). 
This is called the Andreev reflection. In superconductor- 
normal-metal-superconductor junctions, particles (and 
likewise for holes) can be reflected by the two NS in- 
terfaces repeatedly and form bound states in the nor- 
mal region, which are the Andreev bound-states. For 
each transmission channel, there can be a pair of An- 
dreev bound-states which carry currents in opposite di- 
rections. In dc Josephson effects, these bound states are 
responsible for carrying the supercurrent. 4 ' However, it is 
also known that for unequal gap junctions the Andreev 
bound-states can be missing for some ranges of phase 
differenced^ In this case, the supercurrent is then carried 
entirely by the continuum states and the thermal noise 
of the current exhibits dramatically different features^ 
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Previously, ac effects in asymmetric junctions have 
been considered using the scattering matrix approach. 7 
In these works, the authors focused on the character- 
ization of the additional subgap structures due to the 
presence of two superconducting gaps and on the possi- 
ble scheme for the measurement of the phase difference* 7 " 
In this article, we shall instead concentrate on the dy- 
namics of the Andreev bound-states in the ac effects at 
low bias. As pointed out by Averin and Bardas^ in this 
limit the dynamics of Andreev bound-states play a key 
role in the ac effects. Since for unequal gap junctions the 
Andreev bound-states can be missing for some ranges of 
phase difference^ it is therefore of interest to investigate 
its consequence in the ac effect. 

We shall study the ac Josephson effect in unequal- 
gap superconducting quantum point contacts using the 
quasiclassical Green's function method. 9 This approach 
has previously been applied to the study of ac effects 
in symmetric junctions. 10 Here, we shall examine for 
asymmetric junctions the current-voltage characteristics 
at arbitrary transmission coefficients in the point con- 
tact. We will then study in detail the low bias regime 
where inelastic scattering rate is vanishingly small (com- 
pared with the bias voltage and the gap amplitudes) iLll 
Under this situation, similar to the dc effects in unequal- 
gap junctions?*^ the current receives contributions from 
both the Andreev bound-states and the continuum states 
(see Fig. We shall show that in this limit the 

current-voltage characteristics can be understood from 
the current-phase relation provided the different behav- 
iors of the bound states and the continuum states are 
taken into account. For unequal gap configurations, 
where the Andreev bound-states are missing for some 
ranges of phase difference, the continuum states give rise 
to non-zero sine components. Finally we will demon- 
strate that in the low bias limit the temperature depen- 
dence of the dc component is determined from the quasi- 
particle occupations at the smaller gap. In the Appendix 
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FIG. 1: Schematics of the asymmetric junction at low bias. 
The dotted lines at the center of the energy gaps depict the 
chemical potentials of the superconducting banks. The pair 
of discrete states in the normal region are the Andreev bound 
states. The gray areas indicate the continuum states that can 
contribute to the current, for instance, via Andreev reflec- 
tions illustrated by the arrows [with a particle (solid arrow) 
converted into a hole (dashed arrow)]. 



we provide details needed for the calculation. 



II. FORMULATION 

We consider two s-wave superconductors connected by 
a single- mode quantum point contact, which has trans- 
mission probability D. We shall assume that the point 
contact is short compared with the coherence lengths 
of the two superconductors. The order parameters of 
the left and right superconducting electrodes are, respec- 
tively, A; and A r exp(i0), where A^ r are taken to be 
real positive. Without loss of generality, we shall assume 
A; < A r . The junction is brought out of equilibrium by 
connecting the right electrode to a voltage source at fixed 
bias V, while the left electrode is grounded. 

In the quasiclassical Green's function approach to 
superconductivity^ one reduces the Nambu Green's func- 
tion by first separating the fast (relative) and slow (cen- 
ter of mass) degrees of freedom. Since for low energy 
phenomena the relevant time scale is much longer than 
the inverse of the Fermi energy. One can remove the 
irrelevant (fast) degrees of freedom by integrating out 
the magnitude of the relative momentum, retaining only 
the angular information specified by the unit vector p of 
the relative momentum. This leads to the quasiclassical 
Green's function <?(p), which is in general an 8x8 matrix 
in the Keldysh<g)spin(g)particle-hole space (see Appendix). 
Here, for brevity, we have omitted the time and frequency 
variables. For clean superconductors, the Green's func- 
tion satisfies the equation of motion^ 

[ef 3 -A,5]+iv F -V5 = 0, (1) 

where vp is the Fermi velocity, A the off-diagonal self- 



energy (or the pairing function), and 73 the Pauli matrix 
(see Appendix for details). As usual, [a, b] = ab—ba is the 
commutator and all products here involve matrix multi- 
plications and convolutions in energy variables, 9 which 
have been omitted for brevity. Besides the equation 
of motion QJ, the Green's function has to satisfy the 
normalization condition (g) 2 — — 7r 2 l and appropriate 
boundary conditions. 13 Similar to the equilibrium case^ 
the current can be expressed in terms of the difference be- 
tween quasiclassical Green's functions along the incident 
(p) and the reflected (p) directions near the interface 
d = g(p) — <?(p). The quantity d can be solved analyti- 
cally and yieldsi^i^ 



- iD r , / D 



(2) 



The Green's function for the left electrode g l x remains 
its equilibrium form while g r for the right electrode now 
depends on the bias V (see Appendix). 

In the quasiclassical formulation, the current can be 
expressed (we take H = 1, the charge of electron e, and 
the electric current from right to left electrodes positive) 



J(i) 



where TT4 is trace over spin(g>particle-hole space, f z the 
Pauli matrix, and d< = [d K - {d R - d A )]/2 with d*< A > K 
the retarded, advanced, Keldysh components of d, re- 
spectively. Since (d R — d A ) is proportional to the density 
of states, it does not contribute to the total current. One 
can find the explicit form of d K after some algebra and 
then calculate the current lit). These calculations are 
outlined in the Appendix. Choosing the origin of time 
so that the superconducting phase 4> — at t — 0, we 
express the current as a sum over its Fourier components 
in harmonics of the Josephson frequency 2eV 



2ineVt 



(4) 



n— — 00 



where the current components are given by (taking k E 

i) 



In — -^—8n0 



detanh(^) J„(e) (5) 



with i?. 



2 D/tt the normal conductance. Here J n (e) 



is the n-th harmonics of the current density, which has 
a complicated structure given in the Appendix. In the 
tunneling limit (D <C 1), one can easily reproduce pre- 
vious analytical results from these formulas^ Since the 
current I(t) is real, it follows that the current compo- 
nents satisfy I* = Therefore, alternatively one can 
express Eq. Q as 

00 

I(t) = Io + Yj { J « cos(2neVi) + I s n sm{2neVt)} , (6) 
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FIG. 2: The dc (n = 0) current component for asymmetric 
junctions (A;/A r = 0.5) at different transmission coefficients 
(lower to upper curves) D = 0.01, 0.1, 0.4, 0.7, 0.9, 0.99, and 
1 at zero temperature. In all figures we set the parameter for 
inelastic scattering rate r\ = 10 A r . 

where J£ = 2Re{/„} and J£ = 2Im{/„}. In the fol- 
lowing, we shall call the n = 1 components If and If, 
respectively, the cosine and the sine components. 

Applying the formulas obtained above, one can calcu- 
late the current components /„ for general transmission 
coefficient D at arbitrary bias voltages V for any temper- 
ature T. In the following section, we shall demonstrate 
numerical results for the first two current components 
To and h- In particular, we shall study in detail the 
regime where the dynamics of the Andreev bound-states 
is important. This corresponds to the low bias limit for 
which quasiparticle damping is extremely weak. Namely, 
we will be interested in the regime where the inelastic 
scattering rate rj (see Eq. (|A 2(1 in the Appendix) is van- 
ishingly small, so that r\ <C eV <C A; iLiS 



III. CURRENT COMPONENTS 

We shall now study numerical results obtained from 
the formulas derived in the previous section. We will 
first examine the current components at zero tempera- 
ture and then consider the temperature dependence of 
the dc component at the end of this section. In all nu- 
merical results presented below, unless stated otherwise, 
we choose the gap ratio A//A r = 0.5 and take the inelas- 
tic scattering rate rj — 10 _5 A r for the calculations. 

Figure [21 shows the zero temperature results for the 
dc components Iq at different transmission coefficients D 
in the asymmetric junction. As was found previously, 7 
for D = 1 the dc component Jo has a finite interception 
2eA;/-7r at low bias. For bias below (A; + A r ), due to 
multiple Andreev reflections, there are subgap structures 
which are richer than the equal gap case owing to the 
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FIG. 3: The zero temperature results for (a) the real and (b) 
the imaginary parts of the current components 7i for asym- 
metric junctions (A;/A r = 0.5) at different junction trans- 
parencies. 

presence of two gaps. These structures can be classified 
as detailed in Ref. and we shall not repeat it here. 

Figure |3 shows the real and the imaginary parts of the 
current component I\ at zero temperature. The general 
features are very similar to the equal gap case (cf. Fig. 2 
in Ref. |8|). In the present case, however, the current 
components undergo larger oscillations in the subgap re- 
gion than the equal gap results. Moreover, for D = 1 
the imaginary part of I\ (and hence the sine component 
I() can take non-zero values - in contrast to the equal 
gap case, where the sine component vanishes identically. 8 
We shall see below that in the low bias limit this finite 
sine component can be understood by generalizing the 
picture obtained by Averin and Bardas 8 - to the unequal 
gap situation. Namely, we will show that the non-zero 
sine component originates from qualitative changes in the 
current-phase relation when the gap ratio is not one. It 
will be seen that the contribution arises completely from 
the continuum states. 
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FIG. 4: The zero temperature current-phase relations for 
ballistic junctions (D — 1) with gap ratio Ai/A r = 0.5 
(a) in equilibrium (V = 0) and (b) at low bias (V — > 0). 
The total current (full lines) and the current contributions 
from the bound states (dot-dashed lines) and the continuum 
states (dashed lines) are shown. Note that the continuum 
contributions remain the same, whereas there is no branch 
switching for the bound state for V 7^ (see text). The in- 
set in (a) displays the bound-state spectra for the particle 
right-moving (dot-dashed line) and left-moving (dotted line) 
branches. Note that the bound states do not exist for some 
ranges of 



As pointed out by Averin and Bardas for equal gap 
junctions^ at low bias the ac current components can be 
related to the Fourier components of the current-phase 
relation. For equal gap junctions, the supercurrent is car- 
ried entirely by the Andreev bound-states. For unequal 
gap junctions, however, apart from the bound states, the 
continuum states (A/ < |e| < A r ) can also contribute to 
the currenti^& As shown in Fig.0fa)(b), the continuum 
contributions behave the same under zero and finite bias. 
In particular, the continuum contribution is always odd 
with respective to <f> = tt. For the bound-state contribu- 
tions, however, the situation is quite different, as we shall 



FIG. 5: The variation of the first Fourier components at low 
bias with respect to gap ratio for ballistic (D = 1) junctions. 
The lines depict the first Fourier components of the current- 
phase relation, while the symbols are the real (crosses) and 
imaginary (circles) parts of the current component I\. The 
deviations at small gap ratios are due to numerical difficulties 
in low-bias calculations for 1%. Note that for unequal gaps, 
the sine component is no longer zero. 



now explain. 

The two branches of the Andreev bound-states carry 
supercurrent in opposite directions and only one of them 
can be occupied (see inset of Fig. Ela)). In equilibrium, 
since the chemical potential lies in the middle of the en- 
ergy gap (e = 0), it is always the lower branch that is oc- 
cupied. Therefore in this shown in Fig.^a), the 
bound-state contribution always switches branch when 
the superconducting phase <fi goes across 7r. 5 In the pres- 
ence of finite bias, however, since the phase evolves ac- 
cording to the Josephson relation (dc/> /At = 2eV), the 
bound-state current always follows a single branch of the 
bound-state spectra, so that the corresponding current 
contribution remains the same sign during one period. 8 
Figure0Jb) shows the low-bias effective current-phase re- 
lation for ballistic junctions with gap ratio A//A r = 0.5. 
Since there is no branch switching, the bound-state cur- 
rent stays positive and is even with respect to c/> = n. On 
the other hand, due to the finite gap for the continuum 
states, their contributions is the same as the V — case 
and remains odd with respect to cf> = ir. Since sinc/i is an 
odd function with respect to <p — ir, the sine transform of 
the low bias current-phase relation is thus non-zero due 
entirely to the continuum-state contributions. 

To verify this picture, we show in Figure a compar- 
ison between I\ and the first Fourier components of the 
low-bias effective current-phase relation of Fig. 01b) for 
ballistic junctions at different gap ratios. One can ob- 
serve that they match well with each other. The sine 
components become non-zero when the gap ratio is not 
equal to one. As pointed out above, this can be under- 
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FIG. 6: The current densities Jo for asymmetric junctions 
with A;/A r = 0.5. We plot only the region e > since the 
current density can be taken an odd function of energy (the 
even part will not survive the integral in Eq. @). The main 
panel shows the plot for D — 1 at bias eV = 10~ 3 A,-. Inset (a) 
shows current density near the smaller gap for D = 1 (solid 
line) and D — 0.99 (dashed line). Inset (b) shows current 
densities for D — 1 at different bias voltages eV/A r — 10~ 3 
(solid line), 2 x 10~ 3 (dashed line), 5 x 10~ 3 (dotted line), and 
10 x 10" 3 (long dashed line). 

stood from the current-phase relation and is entirely due 
to the continuum states. This result thus confirms that 
the picture obtained in Ref. |8| can be extended to un- 
equal gap junctions. Also, it illustrates the importance 
of continuum-state contributions to the ac current com- 
ponents in unequal gap junctions. 

Finally, we study the temperature dependence of the 
dc component for ballistic junctions in the low bias limit. 
In the equal gap case, it was shown by Averin and 
BardasS that in the low bias limit, the current receives 
major contributions from the neighborhood (~ eV) of 
the gap edges. This is because particles and holes with 
incident energies in these ranges can undergo divergent 
numbers of Andreev reflections (of order A/eV) which 
generate dominant contributions to the current. There- 
fore, it follows that the temperature dependence of the 
current is determined by the occupation factor near the 
gap edge. We shall now generalize this picture to the 
unequal gap configurations. 

To achieve this, we plot in Fig. the current density 
Jo for the dc component at small bias for D = 1. It is 
clearly seen that the current density has a sharp peak 
near the smaller gap A/. Therefore, from Eq. JSJ, one 
concludes easily that the temperature dependence of the 
dc component is determined from the quasiparticlc oc- 
cupation near the smaller gap. This is confirmed by the 
plot shown in Fig. |7| where it is seen that the result fits 
well with the expression (2eA;/7r) tanh(A;/2T). Indeed 
the reason for the sharp peak near the smaller gap is 
very similar to its equal gap counterpart. Namely, at low 



FIG. 7: Temperature dependence of the dc component Io at 
low bias for a ballistic asymmetric junction (A;/A r = 0.5). 
The triangles are results obtained from our expression and 
the solid line is the curve (2eA;/-7r) tanh(A;/2T). Here the 
bias voltage set at eV = 10~ 3 A r . 



bias, excitations near the smaller gap are being injected 
into the point contact and "become" the Andreev bound 
state: they undergo a divergent number of Andreev re- 
flections for V — » 0. For particles and holes with inci- 
dent energies between the smaller and larger gaps, they 
are only partially Andreev reflected and hence contribute 
insignificantly to the current as V — > 0. The presence of 
finite reflection can reduce the current density near the 
smaller gap significantly; this is illustrated in inset (a) 
of Fig. [S] Also, when the bias increases, the peak at 
the smaller gap broadens (see inset (b) of Fig. and 
the current receives contributions from a wider range of 
energies. Eventually, when the bias is large the picture 
described above becomes no longer valid. 

In closing this section, we would like to point out an 
interesting feature which appears in the plots for the cur- 
rent densities. As can be seen clearly from insets (a) and 
(b) of Fig. El there are regular steps (or shoulders) in 
the current density profiles which are separated by 2eV. 
Indeed these steps are consequences of the change in the 
number of Andreev reflections when varying the energy 
of the incident particles. Their occurrence thus further 
supports the picture presented above. 



IV. SUMMARY AND DISCUSSION 

We have studied the ac Josephson effect in asymmet- 
ric superconducting quantum point contacts, where the 
gap amplitudes on the two superconducting electrodes 
are different. Applying the quasiclassical method, we are 
able to derive formulas applicable to general transmission 
coefficient D and arbitrary bias voltages. We calculate 
the IV curves for general D and study in detail the low 
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bias limit, where quasiparticle damping can be ignored. 
We find in this limit for ballistic junctions (D = 1) that 
the sine component becomes non-zero once the two su- 
perconducting gaps become unequal. By comparing the 
results with the Fourier components of the current-phase 
relation, we are able to confirm that the non-zero sine 
component is entirely due to the continuum states. In 
the same limit, we also find that the temperature depen- 
dence of the dc component is determined by the occu- 
pation factor at the smaller gap. These results confirm 
that the picture obtained by Averin and Bardas^ can be 
generalized to the case of asymmetric junctions. 

Finally, to test our results for the non-vanishing sine 
component for unequal gap junctions, it is necessary to 
measure the current component I\ experimentally. As 
suggested by Hurd et al.£ this may be achieved by using a 
proper phase-biasing network. Moreover, in the presence 
of a microwave radiation, it may also have consequences 
on the subharmonic Shapiro steps J^- An interesting issue 
which is not addressed in this article is the change in 
the non-equilibrium noise (shot noise) when the Andreev 
bound-states are missing in asymmetric junctions. This 
will be the subject of our future publication. 
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Here f i are Pauli matrices in the particle-hole space; the 
functions g l and /, are 



v ~ u n 

f^ A (e) = -„-= A ' (A2) 

y/Af - (£ ± IT]) 2 

with r\ a small positive number related to the inelastic 
scattering rate. The Keldysh Green's function is given 
by 

g^gfy-ntf (A3) 
with n l the distribution function 

n l (e,e')=2<K8(e-e')n(e)i, (A 4) 

where n(e) = tanh(e/2T). 

In the presence of bias voltage V, the superconduct- 
ing phase becomes time dependent. The non-equilibrium 
Green's function g can be obtained from the equilibrium 
Green's function g^ through the following transforma- 
tion 

g(t,t>) = S(t)g x (t~t')SHt'), 

= (S(t) o \ (g* \ fsHf) o \ 

I S(t)J \ gi)\ SHf)) ' 

(A 5) 

where 



APPENDIX 

In this appendix we outline the essential elements for 
the calculation of the current components /„. We sup- 
ply the explicit expressions for the Green's functions and 
explain the schemes for the numerical calculation of the 
matrix element d K . Finally, we display the explicit for- 
mula for the current density J n . 



1. The Green's functions 

Since we do not consider magnetic phenomena in the 
present work, the spin degrees of freedom merely intro- 
duce a factor of two (cf. Ref. Il7h . The Green's func- 
tions we shall be dealing with are hence reduced to the 
Keldysh®particle-hole space. Since the superconducting 
electrode on the left remains in equilibrium, the retarded 
(R) and the advanced (A) Green's functions in frequency 
space are (a = i?, A) 

g^ (e, e>) = 2*6 (e - e 1 ) [f.gf® - if,/«(e)] . (A 1) 



S(t) 



e i*W o 
e- i# W 



(A 6) 



with the time dependent phase $(t) = <po/2 + eVt. 

Therefore, for the superconducting electrode on the 
right, expressing f ± = (f x ± if )/2, one can find for a — 
R, A 

^(e,e')= 27r[fJ(e-e')g^(e + r m eV) 

-f + 6(e -e' + 2eV)f*{e + eV)e*° 
+f_6(e -e'- 2eV)f°(e - eV)e^° 1(A 7) 

In this expression, g^' A and f^' A are the same as (|A 2() 
after replacing the subscript I by r there. Similarly, the 
Keldysh Green's function g K can be obtained in the same 
way as (|A 3|l with the replacement of all subscripts I by 
r and using the distribution function 







where n(e) is the same as in QA4[I , 



,(A8) 
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2. The matrix d K 

To calculate the current, we shall need the Keldysh 
component of the quantity d which can be found from 
Eq. © 



where we have denoted H (e,e =F 2eV) = and 
H a (e, e) = 7i^; note that here e = e' + 2meV and that 
we have taken e' as the "origin" to = 0. 

To solve (|Xl3|) . we define the transfer matrices t± so 
that 



-D r 



(A 9) 



for 



where g Q = (g a — g"^ ) and Q 

a = R.,A. As usual, convolution over the energy variables 
is understood in the above formula. To obtain the explicit 
form for d , it is necessary to calculate the quantity Q. 
To proceed, we note first that, using the normalization 
condition (g) 2 = — 7r 2 l, one can find 



2 - D 
2 

2 - D 



1 - 



H- 1 



D 



2n 2 (2-D] 



{&,oo,<t} 



(A 10) 



As usual, {a, b} = db + bd is the anticommutator. Note 
that for brevity here and below we omit the superscripts 
R, A for the retarded and advanced functions. 

Using the explicit expressions for g"^ and g a given 

above, one can express H. in frequency space as 



W(c,e') 



2tt 



S(e - e') H° + S(e - e' - 2eV) H+ 



-S(e -e' + 2eV) H~ 



(All) 



We note that in frequency space 7i possesses the struc- 
ture of a "tight-binding Hamiltonian" Therefore, we 
use the ansatz for its inverse Q 



t^Qm = Qm+1 , TO > ; 

i-Q m =Qm-u m<0. (A 14) 

A 13|) becomes a homogeneous equation, 



For to ^ Eq 

from which one can derive the following recursion rela 
tions for the transfer matrices 



i+ - - 


(^1- 







+ ^m+l^m+l ) 1 TO>0; 



+ ^m-l*ro-l J ' TO < . 

(A 15) 

We solve the transfer matrices numerically by truncat- 
ing the recursion relations (|A 15|) at |to| ~ A r / (eF). This 
means that we are ignoring Andreev reflections when the 
energy is cycled above the gap edges (we choose the larger 
gap A r to improve convergence in the numerical results). 
This is justified since the Andreev reflection coefficients 
decay very quickly above the energy gap. One can thus 
obtain the solution at to = 



Qo = [K 



Hq t Q 



Applying IjA 14fl . we can construct 
to and hence obtain Q R ' A 
results back into IjA 10 
expression for d . One can thus obtain 



(A 16) 



for all values of 



using 
and (|A9j) 



Substituting the 
we can then find the 
Eqs. (@J and 
The explicit form for the n- 



© by making use of © 
th harmonic J n of the current density is given in the 
following section. 



Q(e,e')= J2 ^(e-e' -2meV)Q m (e'). (A 12) 



3. current density 



The equation HQ = 1 immediately leads to For completeness, we provide the explicit expression 

for the current density in this section. We first denote 
Hj,Q m + H^Qm-i + H m Q m +i = S m ol , (A 13) the R 7 A, K components of the commutator 



\$ r ,Si.aXM = ^ ~ + 5(e -e'- 2eF)M>) + S(e - e' + 2eV)M"(e)) 

(A 17) 
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and similarly 

[g R g K + g K g A ) (e, e') = 2tt (s(e - J)N* (e) + S(e - e' - 2eV)N* (e) + 5(e - e' + 2eV)A>* (e)) 



(A 18) 



The explicit forms of the matrix elements M (j ± , N Q ± can be obtained by applying the expressions for and g"^ 

given in Sec. 1 of this Appendix. It is convenient to decompose the Keldysh components M g according to their 
dependence on the distribution function n(e) as 



M o (e) = M o (e) n(e) + M Q (e) n(e + eV) + M Q - (e) n(e - eV) , 
M* ( e ) = M K ° ( e ) n(e) + M*"' (e) n(e - eV) + M*~ 2 (e) u(e - 2eV) , 



+ 



M_ (e) = Af_ (e) n{e) + M_ (e) n(e + eV) + M_ (e) n(e + 2eF) , 



(A 19) 
(A 20) 
(A21) 



and similarly for N . Using IjA 10|) . (|A 12f> . together with (|A 17fl — (|A 21fl in i|A9f) and ©, and then shifting the 
dummy variables so that all occupation factors become n(e), one can obtain J5J and find 



J„(e) 



16?r 2 (2- J D) 



x ^ Tr 4 



f z A(e) 



D 



2n 2 (2-D) 



^2 Tr 4 f B m (e) 



(A 22) 



rn— — oo 



where 



i(e) - M;°(e)Q^(eJ+M^ (e)Q^_ 1 (e n ) + M^(e)Q^ +1 (e„) 



+ M o ( e -eV)Q n (e n -eV) + M_ ( e - eV)Q n+1 (e n - eV) + M_ (e ~ 2eV)Q n+1 (e n ~ 2eV) 
+ M*~\e + eV)Q A n {e n + eV) + M*~\e + e^)Q^_ 1 (e„ + eV) + M*~\e + 2eV)Q*_ 1 (e n + 2eV) 



(A 23) 



and 



B m {e) 



^ KO , , a A 



M(e + 2mey) [iV Q (e)Q n _ m (e n _ rn ) + iV + (e)Q„_ m _ 1 (e„_ m ) + N_ (e)Q n _ m+1 (e n _ m ) 
+ M(e + (2m - l)eV) [n^ (e - e y)Q^ m (e„_ m - eV) + A>* +1 (e - eV) Q^_ m+1 (e„_ m - eV) 
+ A4(e + (2m + l)eF) [jV* _I (e + eV)Q*_ m (e n - m + eV) + N*~\e + ey)Q^_ m _ 1 (e n _ m + eV) 
+ M(e + 2(m + l)eV)N*~ 2 (e + 2eV)^_ m _ 1 (en-m-l) 

+ A4(e + 2(m-l) e y)A>^ +2 (e-2eU)Q^ m+1 (e„_ m+1 ). (A24) 
In the above expressions, we have denoted e„ = e — 2neV and 

A4(e) = M*(e)Q* (em) + M* (e)C_i(e ro ) + ^WCi(^) • ( A2 5) 
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